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' We study the effect on quantum spectra of the existence of small circular disks 

Os ' 

in a billiard system. In the limit where the disk radii vanish there is no effect, 

>>' 

however this limit is approached very slowly so that even very small radii 
have comparatively large effects. We include diffractive orbits which scatter 

^ ' off the small disks in the periodic orbit expansion. This situation is formally 

>'■ 

J^j ' similar to edge diffraction except that the disk radii introduce a length scale 

in the problem such that for wave lengths smaller than the order of the disk 

radius we recover the usual semi-classical approximation; however, for wave 

lengths larger than the order of the disk radius there is a qualitatively different 

behaviour. We test the theory by successfully estimating the positions of 

scattering resonances in geometries consisting of three and four small disks. 
PACS numbers: 0.320.+i, 03.65.Sq 
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The presence of discontinuities in classical Hamiltonian systems implies the necessity of 
a closer study of the quantum mechanics when doing semiclassical periodic orbit theory |I[ 
and has been the theme of numerous recent papers The approach is to study the 

quantum scattering problem near the discontinuity, combine this with classical information 
about classical trajectories away from the discontinuity to find global quantities such as the 
trace of the quantum Green function. In doing so, we maintain the local-global duality 
inherent in periodic orbit theory. In this paper we discuss one class of discontinuity, that 
of small circular scatterers. In the context of the Sinai billiard, the perturbative effect of 



a small disk in the quantum |L0] and classical problems [|TT|JT2[| was studied, but not with 
a scattering interpretation. By small, we mean smaller than the typical wavelength in the 
problem. The opposite limit, of scatterers much larger than a typical wavelength, can be 
evaluated using classical periodic orbits reflecting off the disk plus creeping diffraction to 
account for the discontinuity associated with glancing orbits 0-||,[7|. 

We will analyse billiard systems in two dimensions and therefore seek the Green function 
of the Helmholtz equation 

(V 2 + k 2 ) V = (1) 

with some specified boundary conditions. In the absence of any boundaries, the Green 
function between a source at x' and a receiver at x is 

G f {x,x',k) = - l -H [ +\k\x - x'\) 
i 

e *(fc|x-a!'|-37r/4) _ ^) 



^J%Tik\x — x'\ 

The second line, which is asymptotic in k\x — x'\, will be useful later in the discussion. In 
the exterior of a disk of radius a with Dirichlet boundary conditions, the first line of (0) is 
modified to 



G{xy,k) = -± £ e im ^H^\kx') [ H£\kx) - -^f^HM(kx) ) x> > x 



m=—oo 



G f (x,x',k)+ Z - £ e^-^^^Hlt\kx')H^(kx), (3) 

m=— oo Urn, \KQj) 



where 9 and 9' are the polar angles of points x and x' as measured from the disk centre. 
The first line follows from using Graf's addition formula for Hq + \z) [[13] together with the 
S-matrix of the disk scattering problem. The second line can be seen to equal the first 
by another application of Graf's addition formula and the expansion of J m {z) in terms of 
H^\z). If x > x', one must interchange x with x' and 9 with 9' in the first line of Eq. ([3]) 
but the second line rests unchanged. This has the appealing structure of being the free space 
Green function plus a correction. Whether the correction is small or large depends on the 
geometry of the problem and on the wave number k. 

We now make the assumption that the disk radius a is much smaller than the typical 
distance to points x and x' from the disk. We further assume the semiclassical condition 
kx,kx' 3> 1. Recalling the asymptotic relation H^\z) « exp(— irrnr/2)HQ + \z) (assuming 
\z\ ^> 1 and \z\ > m), we obtain 

£ ^ m{9 - e ' ] ^^H^(kx')HW(kx) ~ 

m=—oo tlm \K(X) 

H^\kx')Ht\kx) £ e^-^-^Pr- (4) 

m=—oo Hm \KQj) 

Although the asymptotic forms used for H^\z) breaks down for m > \z\, this is not 
important in Eq. (Q) since the condition x, x' 3> a implies that the breakdown begins for 
values of m such that the factor J m {ka) / [ka] is already very small. We then conclude 
that in the presence of a small disk centred at position £, the Green function between two 
points separated by an angle <fi as measured from £ is approximately 

G(x, x' , k) w Gf(x, x , k) + d((f>)Gf(x, £, k)Gf(£, x , k) (5) 

where we have defined a diffraction constant 

d(<f>) = -Ai £ e im ^ J r:f q) . (6) 



m=— oo 



Using the relations J- m (z) = (—l) m J m (z) and H^^z) = {—l) m Hl^\z), we can replace 
the exponentials in this expression by cosines thereby displaying the time-reversal property 
d((f>) = d(—(j>). Furthermore, this sum converges since the ratio J m {ka) / H^\ka) decreases 



factorially with m for m > \ka\. In the more general case of a non-circular but small scatterer, 
we must insert the full S matrix into Eq. (|3|). The analysis remains largely unchanged; in 
particular Eq. (^) still applies however the diffraction constant is then a function of both 
angles and not just their difference. We then obtain 

oo 

d(6,6')=2 ]T e i(m'e-me'-( m+ m')n/2) Tmm ^ ( ? ) 
m,m'=— oo 

where we have defined the S matrix through S mm i = 8 mm i — iT mm i . 

The Green function (gD has a direct contribution as if there were no disk plus a contri- 
bution in the form of a product of Green functions which arises from scattering off the disk. 
This is the same structure which exists in the presence of vertices [|14]] where we obtain the 
diffraction constant from the solution of the wedge scattering problem solved by Sommer- 



feld ||15|| . Despite the similar form, there are two aspects of the problems which are quite 
different. The small disk diffractor has no internal orientation but does have an internal 
length scale, a. In contrast, a wedge has an internal orientation, as given by the direction 
of its normal, but has no internal length scale. The existence of an orientation means that 
there are choices of incoming and outgoing angles for which the vertex diffraction constant 
diverges whereas this never happens for the disk. On the other hand, the lack of an internal 
length scale means that the vertex diffraction constant is independent of k whereas for the 
disk it is clearly /c-dependent. The systems do share the property that we can trivially 
extend them to include problems with a potential V(x) ||. We assume that the potential 
does not change much in a wavelength and thereby compute the diffraction constant with 



k = ^2m(E — V(£))/h, the wave number at the disk. Eqs. (||) and (0) then apply where 
Gf is the Van-Vleck approximation to the Green function for that potential in the absence 
of a disk. 

So far we have made no assumption on the value of ka. For b > 1 we recover the 
expected geometrical structure, as we discuss below. There is qualitatively different be- 
haviour for ka <C 1 and a cross-over for ka ~ 1. In the limit ka — > we note that 
J m (ka) / H^' (ka) ~ in(ka/2) 2m jm\{m — 1)! for m^O and only the m = term contributes 



significantly. We call this the s-wave limit. From the approximations 

J Q (ka) » 1 Y Q (ka) » 2 - (tog + 7e ) (8) 

valid for small ka (where 7e = 0.577. . . is Euler's constant), we derive the s-wave approxi- 
mation 

d ~ : tjx -> (9) 

which is independent of scattering angle 0. As fca — > 0, the denominator of Eq. (g) grows 
logarithmically so that the diffraction constant goes to zero and the disk has no effect, which 
is reasonable. However, this happens very slowly so that even for very small values of ka 
there is still an appreciable effect, as we will demonstrate. 

It might seem surprising that the diffraction constant vanishes as a — > since we are de- 
manding that the wave function vanish at a point, and it might be thought that this should 
have some effect. That this is not so can be understood with the example of an annulus 
in which the central disk is very small. Although the eigenf unction does indeed vanish on 
the disk, it increases very rapidly so that within a small distance the eigenfunction is indis- 
tinguishable from one in which there were no central disk. In this sense, the wavefunctions 
(and eigenvalues) are virtually indistinguishable from those corresponding to the disk-free 
system. In Ref. ||16|| , the author argues that disks of zero radius continue to have an effect, 



but the system he was considering was equivalent to an infinitely thin line charge in an 
electro- magnetic wave guide (see also Ref. [0 for the three dimensional generalisation.) In 



his language, our disk is uncharged and there is no contradiction between his conclusions 



and ours. The difference is in the order one takes the limits a — > and k — > oo. In Ref. [[16 
one starts with the first limit (while maintaining a finite interaction) whereas we consider 
the second limit while holding a small but fixed. As a result, the problems are quite distinct. 
For example, in the short wavelength limit the disks considered here will start having a large, 
classical effect (i.e. when ka is of order 1) whereas in the system mentioned above the effect 
of the scatterer vanishes for large k HI"8f . 



In a system composed of several scattering centres, the Green function between points 
labelled x' and x will receive contributions from paths which diffract several times. The 
contribution of one such path, labelled p, is a simple generalisation of the second term of 
(H) and equals 

G p {x,x\k)nG f {x,£ n ,k)d n ^f[G f (^ (10) 

The quantities £j indicate the scattering points - in this case the centres of the disks with the 
subscripts j indicating the order in which they are encountered for that path. The path is 
composed of n scatterers with a diffraction constant dj associated with each one. The details 
of the trace integral are worked out in Refs. nere we simply sketch the derivation. 

We stress that all of these calculations are to leading order in H (or 1/k in this context.) The 
criterion of stationary phase selects periodic orbits which are everywhere classical except at 
the singularities where they diffract by an arbitrary angle. A periodic orbit is composed of 
n classical segments connected by n diffractions. 

If we identify x with x' in ( |I0"D we see that it is a closed cycle of Green functions in which 
the segment between the scatterers £ n and £1 is "cut" by the point x which lies between them. 
In general, the trace integral associated with a periodic orbit labelled 7 will be evaluated 
by integrating over all choices of x. This means that we must allow x to cut open the cycle 
of Green's functions between any two consecutive scatterers - each possibility physically 
corresponds to x lying between that pair of scatterers. For this purpose, we define a parallel 
coordinate z which runs along the periodic orbit from £ n to £1 to £2 etc. until it returns to 
£ n . When z is between & and £i+x, the point x is between these two scatterers (where we 
identify n + 1 with 1.) At each point along the orbit, we define a transverse coordinate y 
so that x is parameterised by the pair (y,z), as for geometric orbits The integrand of 
the trace integral is similar to flip]) but where we cut the cycle of Green's functions between 
£j and (as governed by z) rather than between £ n and £1 so that the trace integral 
associated with 7 is 
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Due to stationary phase, only those points close to the orbit contribute significantly so that 
to leading order, all the diffraction constants dj are independent of x and can be considered 
invariant properties of the periodic orbit. Therefore, the only x dependence is in the two 
Green functions connecting x to the adjacent scatterers; all other factors in the integrand are 
constant. For each value of z, we evaluate the transverse integral / dy by stationary phase; 
this yields a factor which is independent of z and is proportional to the Green function 
between the two adjacent scatterers. Combining this with all the constant factors in the 
integrand gives the product of the closed cycle of Green functions and diffraction constants 
from one scatterer to the next, IT/=i Gf(£j+ii k)dj. We get this same invariant factor after 
doing the y integral at any point z along any segment of the periodic orbit so there is no 
explicit z dependence remaining. The integral parallel to the orbit § dz = L 7 is then simply 
the length of the periodic orbit, as also happens for geometric orbits. The final result is 

g^^-i^luGfi^^d} 

W -^15-7=4 exp{.(A;L 7 - 3n 7 7r/4)}. (12) 

We have made use of (§) where Lj = — This expression involves one fewer Green 
function than fllQ|) . (The additional Green function contributes a constant proportional to 
e -i37r/4 1 ^jj, whjch combines with a factor proportional to e l7r//4 / y/k from the stationary phase 
integral to give the prefactor of —i/2k.) Eq. flT2|) is the contribution of a single diffractive 
periodic orbit, in general we must sum over all such orbits as well as over all purely geometric 
orbits to get the total trace. For this reason, we have introduced the subscript 7 on the 
index n 7 in the above equation. 

Following Refs. [§,f|,§| we write down the semiclassical diffractive zeta function [JBJ whose 
zeros approximate the exact quantum resonances, 

7 



Cdif = IK 1 - as) 

7 

where 

{ n d ) 

t 7 = < II 3 ( exp{i(fcL 7 - 3n 7 7r/4)} (14) 
This results follows from the semiclassical approximation 

§ - ^ < 15 > 



so that the sum over all diffractive orbits in Eq. (|12|) is the logarithmic derivative of the 
zeta function (0). (We take the derivative with respect to k 2 since the trace of the Green 
function (|i~2|) is properly thought of as a function of k 2 and not of k.) The product is over just 
the primitive orbits; their repeats have already been summed. In a system with coexisting 
geometric and diffractive orbits, we need to multiply the corresponding zeta functions ||[|] . 
The result is a purely formal product which must be regulated differently for scattering 
|2(1|^22| and bound [23[] problems so that its zeros are the semiclassical eigenvalues of the 
full problem and not the zeros of the individual terms in the product. The diffractive zeta 



function ( [[3D involves no additional product as happens for geometric orbits [p^j , resulting 
in there being only leading resonances in scattering calculations [|5]|§. 

We now specialise the discussion to scattering geometries featuring three and four small 
disks arranged symmetrically in the plane |25H32f. We first discuss the three disk problem 



as shown in Fig. |T]a (we exaggerate the size of the disks to make the discussion clearer.) 
Starting from one of the disks, there are two distinct processes. We can go to one of the 
other two disks and either scatter back to the original disk or scatter on to the third disk. We 
assign these two processes the symbols and 1 respectively [p9f| . The weights are t = d(0)u 
and t\ = d(n/3)u where the factor 

u = -fTTB e M*(kR - 3vr/4)} (16) 
y/SnkR 

is common to both orbits and R 2> a is the inter-disk spacing. 



S 



Notice that there is a C^ v symmetry to this problem f27H30|l consisting of the identity, 
rotations by ±27r/3 and reflections through the three symmetry axes. This group has three 
irreducible representations which are called Ax, A2 and E. We can make use of this symme- 
try by considering dynamics in the fundamental domain ||33|| , which in this case is a wedge 
consisting of one sixth of the plane as indicated in Fig. |T]a. One does this by following a 
trajectory and using the symmetry operations to map the trajectory back into the funda- 
mental domain whenever it crosses a boundary. In the fundamental domain of the three 
disk problem, there is only one half disk. A trajectory can leave this half disk in only one 
direction, which is labelled A. Upon encountering the border of the fundamental domain, a 
reflection operation is applied so that the trajectory returns to the disk, where it has two 
choices. It can either diffract back onto A or it can diffract into the direction A'. In the 
second case, we apply a reflection operator again to map this back onto A. These two possi- 
bilities are both diffractive periodic orbits of the fundamental domain and have the weights 
t and £1 discussed above. Each orbit has an additional group theoretic weight given by the 
characters (in the representation being considered) of the group operations needed to keep 
the orbit in the fundamental domain |30| . 

In general there are longer periodic orbits as labelled by whether they back scatter or 
forward scatter at each encounter with the disk. These can be then labelled by a binary 
sequence of O's and l's. However, there is a multiplicative property to the weights such 
that the weight of any long orbit is equal to the product of the weights of shorter cycles. 
For example tool = ^0^1 si nce they both equal d^diu 3 . This property means that we can 
represent the zeta function as being the determinant of a Markov graph p4 |, which is drawn 
in Fig. [L]b. The single node in the graph, A, is connected to itself by the processes and 1 
described above. 

All the characters of the totally symmetric representation Ai are unity, which simplifies 
its discussion. To find its zeta function, we simply read off from the Markov graph all 
nonintersecting closed loops. In this case there are only two and we get the simple result 
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Gl = l-t -h. (17) 
This formula agrees with the result found in Refs. p9|j30| for the special case where all the 



higher order "curvature corrections" 20-22 vanish identically. This vanishing is simply 



a result of the fact that we have a one node graph so we only need consider weights of 



topological length one. Armed with this rule, we can then read off from Refs. |29|,|3(J the 
zeta functions of the other two representations. These are 

Cm = i + *o-*i Ci 1 = i + *i + *?-*o- (18) 

We could have ignored the symmetry decomposition and simply drawn the six node Markov 
graph of the full problem as shown in Fig. |2|. The vastly increased number of closed loops 
in comparison to Fig. [3] underlines the advantage of using the symmetry reduction. The 
rule for finding the zeta function is to find all non-intersecting closed loops and products of 
non-intersecting closed loops. A product of n non-intersecting closed loops has a relative 
sign (— l) n . Carefully enumerating all such loops of the full graph, we find its zeta function 
to be 

C 1 = 1 - St 2 - 2t\ + 3t 4 - 3t 2 tl -4 + 4- 344 + 3to*i- ( 19 ) 

This equals the product CXl Ca 2 Ci 2 °f the symmetry decomposed zeta functions above. In 
addition to the additional complexity of its Markov graph and zeta function, the full zeta 
function has the further disadvantage that we do not know to which symmetry class one of 
its zeros belongs. However, this exercise is useful in verifying that our use of the results of 



Ref. [|9],[3g is well founded. 

The exact resonances of this geometry can be found numerically by finding the zeros of 
the determinant of a matrix. This matrix is [ p8[ 



M nm — 5 nm + A nm (20) 

where for the A\ resonances 
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_ J n {ka) 



cos (|(5n - m)) H { +J m (kR) + (-) m cos (|(5n + m)) H { n + + UkR) 



21) 



#i +) (A;a) 

Expressing detM in a cumulant expansion @JS^|, valid because A is trace-class yields 



detM = 1 + trA - ~ (trA 2 - (trA) 2 ) + • • • (22) 
where from Eq . (|2l|) one obtains 

^ = £ ^^(cos(^)^ +) (^)+^: ) (^)). (23) 



m=— oo 



We now impose the same constraints as before, namely kR 3> 1 and R ^> a so that we can 
replace H^(kR) by cos(mir)HQ + \kR) and 

tiA^H { + \kR) £ J ^ fca) (cos(27rm/3) + cos(mTr)) . (24) 



m=—oo 



Then using Eqs.(^) and fll6|) and the asymptotic form of Ho~\kR), we find 

trA « — it (d(0) + d(7r/3)) = -(f + h). (25) 

We see that truncating the cumulant expansion at the term linear in A and invoking the 
relevant approximations gives the same equation for detM = as we earlier derived for 
= 0. This is reassuring since it means that we understand the error caused by replacing 
H^\z) by exp (— irmc /2)Hq + \z) in Eqs.(£D and (|24]) ; it is the same as neglecting higher 
order terms in the cumulant expansion. As shown in Refs. P,|35|, this is equivalent to 



neglecting higher order curvature corrections in the cycle expansion. The fact that the 
semiclassical approximation can be made on the level of the traces of the scattering kernel 



trA n (|21|) which result from the defining cumulant expansion ( pffi) provides an alternate 
method to arrive directly at the zeta function Cjjff- This method does not require closed 
expressions for the trace of the Green function, does not invoke the semiclassical relation 
(|15|) and, most importantly, appears in a curvature-regulated form In particular, 

one can use this to read the weights t 1 (|TID directly from the Green's function product 



n?=i k)dj, which is just the closed path equivalent of the open path Green's 



function (II 
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The identification between the quantisation conditions Ca.i = and detM = tells us 
something else. In Ref. 0] it is shown that one can extract the contribution of geometric 
orbits and diffractive creeping orbits from tiA by invoking Watson contour integration to 
replace the sum of Eq. fl23|) . This means that the diffraction constant contains information 
about periodic orbits and creeping. Therefore, even in the limit ka ^> 1, the formalism 
described here still applies, the price being the necessity to include many terms in calculating 
the diffraction constant @. We therefore have a uniform picture. For large values of ka, 
one invokes geometric and creeping orbits but for intermediate and small values one invokes 
the small disk scattering theory elucidated here. These are guaranteed to match smoothly. 
Although this was shown explicitly only for two and three disk systems, the same will hold 
for any number of disks in any geometrical arrangement. 

We show the exact and semiclassical results in Fig. ^|for the A\ and E resonances together 
with the approximations using geometric orbits. The resonances are shown in the complex k 
plane and are measured in units of 1/R. In Fig. |3|a we show the results for the A\ resonances 
for R/a = 60 so that the cross-over condition Realjfca} = 1 corresponds to R,eal{kR} = 60. 
The minimum, which is developing at the right of the figure, has a geometrical interpretation 
in terms of interference between the two shortest geometrical orbits in the fundamental 
domain, t and t\ [JZ7|,fZ[J. As promised, the diffractive picture captures this behaviour. For 
the highest values of k, we used 70 partial waves in the calculation of the diffraction constant 
(0). If we held the number of partial waves fixed, the calculation would start to fail for larger 



values of \kR\. We also include the results from the theory of geometrical orbits |29] for 
comparison. The new regime is at the left of the figure where Realjfca} <C 1. There it can 
be seen that the widths of the resonances increase logarithmically with kR, a result which 
we generically expect for diffraction |5]||. In those references it is shown that the width of 
the first resonance scales as log(ci) and since d scales logarithmically with a, we find that 
the width of the first resonance (when measured in units of 1/-R) scales as log(log(i?/a)), 
as opposed to the \og(R/a) behaviour predicted by geometric orbits H. This means that in 
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the diffractive case, the resonances are observable even for extremely large values of R/a. 

In Fig. |3|b we show the results for the A\ and E resonances for Rj a = 600. The agreement 
conforms to the discussion of the top panel, however the increased value of R/a means that 
none of the resonances shown are in the geometrical regime Real{/ca} > 1 so there is no 
strong interference between the weights t Q and t\. The A% resonances have smaller widths 
because, to linear order, their zeta function (|i~7f) involves two weights, which are in phase, 
while the zeta function for the E resonances fll8|) involves only one weight. A disadvantage 
of the three disk problem for this study is that in the diffractive regime to C 1 the s-wave 
term dominates so that the two diffraction constants d(0) and d(ir/3) are very nearly equal 
and so too are the weights to an d t%. The result on the spectrum is approximately the same 
as if there were just one weight, a situation which is known to lead to rather uninteresting 
spectra For this reason we were led to study the four disk problem which we discuss 

next. 

The four disk problem shown in Fig. ^a has more structure than the three disk because 
there are two distinct lengths in the problem; in addition to the side length R, there is the 
diagonal length y/2R. Accordingly, we define the factor 

1 exp{t(y^2fci?-37r/4)}, (26) 
8V2nkR 

in analogy to u. In addition, there are now six distinct processes. Starting at any disk we can 
go to one of the two near disks and either diffract back with do = d(0), diffract to the next 
disk with d\ = d{it/2) or diffract diagonally with d 2 = d(ir/A). Additionally we can head 
diagonally across and diffract back with d^ = d(0) or to either one of the other two disks, 
again with d^ = d(n/4). This problem has C± v symmetry which has four one dimensional 
representations labelled A±, A 2 , B\ and B 2 and one two dimensional representation labelled 
E [BO]. This system has previously been studied semiclassically using periodic geometric 



orbits |32|] 



As before, we want to find the Markov graph of the problem for which we study the 
dynamics in the fundamental domain which is one eighth of the full plane and is shown in 
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Fig. |]a. Starting at the half disk, we can go in one of two directions, which we call A and B. 
We want to find all paths which start and end at either A or B. From A we first reflect off 
the diagonal wall and upon returning either diffract back which we call 0, diffract to A' and 
then reflect onto A which we call 1 or diffract to B which we call 2. From B we first travel 
to the centre and on returning either diffract back which we call 3, diffract to A which we 
call 4 or diffract to A' and reflect to A which we call 4. This is shown diagrammatically as a 
Markov graph in Fig. f§b. Note that process 3 is a boundary orbit which lies on a symmetry 
axis and can be shown to contribute only to the spectra of representations which are not 



odd with respect to reflections through that axis |36 . 

The weights corresponding to each process involve one geometric arc and one diffraction 
so we find 

to = dou t\ = d\u t 2 = d 2 u t 3 = d 3 v t± = t| = d±v. (27) 

In general each one of these also has a group theoretic factor depending on the group 
representation being considered. Again, we start with the symmetric A\ representation for 
which all the characters equal one. Enumerating all closed loops and products of closed 



loops on the graph, we read off the zeta function |L9 



as 



Cm = 1 - to - ti - * 3 - (2t 2 t 4 - t t 3 - tit 3 ) (28) 

where we have used the equality between t 4 and t^. This result involves cycles of topological 
lengths one and two. We now have contributions of length two since the graph has two 
nodes, however cycles of length three and higher are absent in Eq.(^). We again note 
that this is the same expression as the cycle expansion of the 4-disk problem discussed in 



Ref. |30| where we use use toi = t ti and additionally invoke the identification between 
{t 0} ti,t 3 , t 2 U, t 2 t 4} tit 2 U, t t 2 t 4 } in our notation and {t , tx,t 2 , t 02 , t l2} t n2 , t 002 } in theirs. As 
before, we can use this fact to read off the zeta functions of the other representations from 



Ref. [30 



Cb 2 — 1 + ^0 + ^1 — ^3 + (2^4 — ^0^3 — 
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= 1 + 1 - h CbI = 1 - *o + h 

G 1 = 1 + h + {t\ - to) + (2*0*2*4 - tlh - 2*i* 2 * 4 + t%). (29) 

We are primarily interested in the region ka < 1 for which the diffraction constants are 
almost equal (ie the s-wave limit) and we see that the A 2 and B\ representations have almost 
total cancellation and therefore their resonances are comparatively deep in the complex 
k plane. These are the representations which are odd with respect to reflections across 
the diagonals of the square so process 3 does not contribute. Instead we concentrate on 
the representations A±, B 2 and E. In Fig. [| we plot the exact positions of these three 
representations found using the algorithms of Refs. p5| , p7| together with the semiclassical 
approximations from Eqs. ( p8|) and (p9|) for R/a = 600. In all cases, the semiclassical 
predictions from the zeta functions work well although it is interesting to note that there is 
a noticeable deterioration of the quality for the resonances with large imaginary part. The 
irreducible representations Ai and B 2 have richer spectra due to the interferences among 
the three basic weights. The E resonances are given by a zeta function which is dominated 
by the weight t% and thereby shows the characteristic logarithmic behaviour discussed above 
and observed in Fig. [5]a. For larger values of k, the quadratic terms of Eq. fl2"9"|) become 
important leading to more structure in the E spectrum. This structure will eventually 
develop into the rich spectrum of scattered resonances predicted by the geometrical orbits. 

In this problem, we have altogether defined 5 weights - however since t 2 and *4 always 
occur as a product in the zeta functions, it is more precise to say there are 4 independent 
quantities. On the other hand, it is known that the geometric orbits can be labelled with 
just three symbols |30|. (In contrast, for the three disk problem we had only two weights, 



in agreement with the two symbols needed to label geometric periodic orbits.) This is 
reminiscent of the approximations of the transfer operator based on so-called T matrices 
|38| which lead to transcendantal quantisation equations like (p9|) but in terms of classical 



trajectories. Increasing the dimension of the T matrices induces more complicated equations 
in terms of which the quantisation is more exact. That approximation is based on assuming 
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certain matrix elements (weights in our language) are approximately multiplicatively related 
and so drop out of the equations - as happens exactly for diffraction. The structural similarity 
between these results is presumably based on the underlying structure of finite Markov 
graphs which are used by us and are implicit in the work of ||38|| . 

In conclusion, we have discussed a form of discontinuity which is amenable to discussion 
in terms of diffraction, that of small disks. Since the effect of a disk vanishes as the disk 
radius goes to zero, we must consider disks of some fixed size. Doing so introduces a length 
scale in the problem such that if ka ^> 1 one can use standard geometrical orbits. However 
in the domain ka <C 1 a qualitatively new physical picture is necessary. The formalism 
we discuss here incorporates both limits but at the price of having to include many partial 
waves when fca>l. We have tested this theory in systems consisting of three and four disks 
arranged symmetrically on the plane. The formalism of Markov graphs and zeta functions 
applies equally well to any system in which there exist objects which can be approximated as 



point singularities, including point scatterers mentioned above and Aharonov-Bohm flux 
lines P5 |. These systems allow a finite approximation based on zeta functions to give their 
scattering resonances and as such are formally useful in testing the formalism. However, 
the arguments developed here apply equally well to bound systems. Putting a small disk or 
other singularity inside a billiard introduces diffractive paths which appear in the Fourier 
transform of the spectrum [[EJ in a characteristic way, just as with edge diffraction HH]. 
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FIGURES 

FIG. 1. Left: The configuration space of the three disk problem with the fundamental domain 
indicated at the top right. The arrow A indicates the unique direction in which a trajectory can 
leave the disk and ultimately return and the arrow A' is its mirror image. Right: The corresponding 
Markov graph with a single node A and the two processes which connect it to itself. 

FIG. 2. The same as Fig. 1 but without using the symmetry decomposition. There are now six 
possible directions and consequently a six node graph. The six short lines correspond to weights 
to and the six long lines correspond to weights t\. 

FIG. 3. Top: The A\ three disk resonances for R/a = 60 plotted in the complex k plane. The 
exact resonances are represented as open circles, the semiclassical diffractive predictions as vertical 
crosses and the semiclassical geometric predictions as diagonal crosses. Bottom: The upper set 
of points are the A\ resonances for R/a = 600 while the lower set are the E resonances. Same 
symbol convention but we do not include the geometrical orbit predictions. In both cases the wave 
numbers are measured in units of l/R. 

FIG. 4. Left: The configuration space of the four disk problem with the fundamental domain 
indicated. The two available directions are A and B and A' is the mirror image of A. Right: The 
corresponding Markov graph with two nodes A and B and all the interconnecting processes. 

FIG. 5. Top: The A x resonances of the four disk problem for R/a = 600. Middle: The B 2 
resonances. Bottom: The E resonances. (Same symbol and unit convention as in Fig. 3.) 
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Fig. la 



Fig. lb 




Fig. 2b 



